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We investigate matter-enhanced Mikheyev-Smirnov-Wolfenstein (MSW) active-sterile neutrino 
conversion in the ^ Vs channel in the collapse of the iron core of a pre-supernova star. For values 
of sterile neutrino rest mass rus and vacuum mixing angle 6 (specifically, 0.5 keV < ras < lOkeV 
and sin^ 26 > 5 x 10~^^) which include those required for viable sterile neutrino dark matter, our 
one-zone in-fall phase collapse calculations show a significant reduction in core lepton fraction. This 
would result in a smaller homologous core and therefore a smaller initial shock energy, disfavoring 
successful shock re-heating and the prospects for an explosion. However, these calculations also 
suggest that the MSW resonance energy can exhibit a minimum located between the center and 
surface of the core. In turn, this suggests a post-core-bounce mechanism to enhance neutrino 
transport and neutrino luminosities at the core surface and thereby augment shock re-heating: (1) 
scattering-induced or coherent MSW Ve Vs conversion occurs deep in the core, at the first MSW 
resonance, where I'e energies are large (~ 150 MeV); (2) the high energy Vs stream outward at near 
light speed; (3) they deposit their energy when they encounter the second MSW resonance Vs 
just below the proto-neutron star surface. 

PACS numbers: 14.60.Pq,95.35.+d,97.60.Bw,98.80.-k 



I. INTRODUCTION 

The nature of the dark matter and the explosion mech- 
anism for core collapse supernovae are among the most 
vexing outstanding problems in modern astrophysics. In 
this paper we point out a curious link between these 
problems. We show how sterile neutrinos with rest 
mass and vacuum mixing parameters in a range which 
makes them viable cold and warm dark matter candi- 
dates HIlllllllllllllliBo] also can alter energy 
and lepton number transport physics in collapsing stel- 
lar cores, possibly completely altering the standard core 
collapse/explosion model. 

In broad brush, the electron fraction Yg is lowered, 
lowering the pressure and, hence, the homologous core 
mass. This augers against an explosion because the ini- 
tial bounce shock energy is reduced, leading to the shock 
stalling closer to the neutron star where shock re-heating 
by neutrino energy deposition will be less effective. How- 
ever, we also point out here that alterations in the rate 
and nature of neutrino energy transport engendered by 
active-sterile neutrino flavor conversion could go in the 
direction of helping the explosion by increasing neutrino 
energy luminosities at a critical epoch when the fate of 
the supernova shock is determined. 

When large magnetic fields are present, active-sterile 
neutrino flavor transformation involving sterile neutrinos 
with rest mass ^ keV can facilitate the generation 
of asymmetric neutrino emission from the neutron star 
surface. This, in turn, can generate large "kicks," which 
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may be required to explain the space motions of at least 
some pulsars [nl [12. Recent work suggest that these 
kicks can augment the hydrodynamic transport of neu- 
trino energy to the base of the stalled supernova shock, 
thereby increasing shock re-heating and helping to gen- 
erate a supernova explosion|l3j|. 

The prospect of sterile neutrinos aiding the supernova 
explosion is a concept that challenges the conventional 
wisdom that neutrino energy diverted into sterile states 
is energy unavailable to aid shock re-heating. However, 
it must be kept in mind that the energy (kinetic plus 
optical) in the explosion, ~ 10^^ ergs, is only ^ 10% of 
the neutrino sea energy immediately after core bounce, 
and only a miniscule ~ 1% of the ^ 10^^ ergs in the 
neutrino sea a few seconds after core bounce (tpb > 2 s). 

Perhaps energy is not the problem; but the transport 
of neutrino energy well may be the crux problem. Here 
we show that if neutrinos spend some of their time in 
sterile states, their transport mean free paths can be re- 
normalized upward, thereby increasing the neutrino lu- 
minosity at the neutrino sphere. A modest increase in 
neutrino luminosity during the epoch when the stalled 
bounce shock is being re-heated (ipb < 1 s) may be the 
difference between an explosion and a dud. Though the 
existence of light sterile neutrinos may be unlikely, the 
tremendous leverage they would have in altering core lep- 
ton numbers and transport properties suggests that we 
set aside theoretical prejudice for the purpose of explo- 
ration. 

Active-active channel neutrino flavor mixing, which we 
know exists, may also affect supernova physics. Indeed, 
t^e ~* i^fi.T flavor transformation in the in-fall phase of 
stellar collapse was investigated 14] shortly after the 
discovery of the Mikheyev-Smirnov-Wolfenstein (MSW) 
llSl lid mechanism for in-medium enhancement of neu- 
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trino flavor conversion. However, because of the large 
lepton degeneracy and low entropy in the collapsing core, 
only un-physically large active neutrino masses can un- 
dergo significant matter-enhanced flavor transformation 
there. 

The mass-squared differences, and two of the four 
vacuum mixing parameters, relevant for the 3-active- 
neutrino (i/p. i/| , ,yT- ) case now are known accurately from 
experiment |l7l Il8| . The mass-squared differences are: 
the "solar neutrino scale" 5m? w 8 x 10~^eV^; and the 
"atmospheric scale" 5m? w 3x 10^^ eV^. These are small 
on the scale of neutrino mass-squared difference neces- 
sary to influence core in- fall epoch physics 0| . 

However, the anomalous result of the LSND experi- 
ment ll9|,|20|,|2l|, if interpreted as arising from vacuum 
neutrino oscillations, would seem to suggest the existence 
of a fourth neutrino which, on account of the Z'^-width 
limit, must have interactions considerabley weaker than 
those of the weak interaction. This is then a so-called 
"sterile neutrino," perhaps an SU(2) Standard Model sin- 
glet. If true, this is an astounding result which calls into 
question the "See-Saw" model for how very massive right- 
handed neutrino states force neutrino states closely asso- 
ciated with active flavors to be light. The mini-BooNE 
experiment [22] is currently probing the same parameter 
range as the LSND experiment (and more). 

Regardless of the outcome of this experiment, sterile 
neutrinos with rest masses, and vacuum mixing angles 
with active species, in the ranges 0.5 keV < mg < 15keV 
and sin 20 > 10"^^ respectively, remain a possibility. 
Currently the only constraints on this speculative sector 
of particle physics come from X-ray astronomy 23, 24] . 
The existence of neutrinos with these mass/mixing pa- 
rameters could subtly or, more likely, drastically alter 
supernova core and shock re-heating physics as well as 
the supernova neutrino signal. Lighter sterile neutrinos 
have been considered in other aspects of the supernova 
environment miHlllIll El. The central point of this 
paper is that to assess the effects of these sterile neutrinos 
on supernovae we must include them starting from the 
outset of core collapse, as many "down-stream" aspects 
of composition, equation of state, and energy/entropy 
transport which affect the explosion mechanism, nucle- 
osynthesis, and expected neutrino signal could be af- 
fected. 

In Section II we discuss active-sterile neutrino flavor 
conversion through MSW resonances in the collapsing 
core. We examine the prospects for energy and lepton 
number transport enhancement in Section III. We give 
a discussion and conclusions in Section IV. Appendix A 
describes the liquid drop model equation of state and the 
lepton number re-distribution schemes employed in our 
one-zone collapse code. 



II. RESONANT ACTIVE-STERILE NEUTRINO 
FLAVOR CONVERSION IN STELLAR 
COLLAPSE 

A. The Standard In-Fall Epoch Collapse Scenario 

The weak interaction in general affects nearly every 
aspect of the core collapse supernova phenomenon. Neu- 
trino emission, scattering, and absorption signiflcantly 
influence the core collapse epoch and the prospects for 
obtaining a viable supernova explosion. It is therefore no 
surprise that massive sterile neutrinos could grossly alter 
the standard model for supernova core collapse. 

Stars with initial masses > 10 Mq evolve through a 
succession of nuclear burning phases, dumping entropy 
along the way through neutrino emission, eventually 
forming a Chandrasekhar mass (~ 1.4 Mq) core sup- 
ported by relativistically degenerate electrons, possessing 
a low entropy s « 1 (in units of Boltzmann's constant per 
baryon), and composed of iron-peak nuclei in nuclear sta- 
tistical equilibrium (NSE). This core will go dynamically 
unstable and will collapse to nuclear density on a time 
scale ~ 1 s. 

Low entropy, electron degenerate conditions determine 
the course of the weakly interacting sector during the 
collapse. Electron capture, which is the forward reaction 
in 

e~+p^n + i/e, (1) 
proceeds rapidly because the electron Fermi energy is 

1/3 

large, /Lte « ll.lMeV{pioYe) ' , where Ye is the number 
of electrons per baryon and pio is the density in units 
of 10^°gcm~^. Because the entropy is low, the proton 
targets for electron capture reside principally in heavy 
nuclei. At first the ^'e's produced by electron capture es- 
cape and carry away entropy. However, electron capture 
can leave daughter nuclei in excited states with excitation 
energies well above the temperature T, thereby increas- 
ing the entropy per baryon. On balance, the entropy 
will rise modestly, but will always remain low, s '--^ 1, 
while the temperature will be in the range T « 1 MeV to 
2 MeV (Here we will measure entropy s in units 

of Boltzmann's constant per baryon.) 

As the collapse proceeds and the density increases, Yg 
will decrease and the mean nuclear mass A will rise. 
Coherent neutral current neutrino scattering on nuclei 
has a cross section that scales like A^, so eventually the 
I'e's produced by electron capture will be trapped - their 
mean free paths will drop below the core size or, equiv- 
alently, the neutrino diffusion time will exceed the col- 
lapse time scale. Neutrino coherent neutral current scat- 
tering on heavy nuclei is essentially conservative. As a 
result, neutrino-electron and neutrino-neutrino scatter- 
ing as well as de-excitation of hot nuclei into neutrino- 
antineutrino pairs must intervene to cause the Vf, dis- 
tribution function to evolve into a thermal, Fermi-Dirac 
form. 
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Eventually all lepton species will come into thermal 
and chemical equilibrium with the nuclear component. 
This is beta equilibrium. When beta equilibrium obtains, 
the rates of the forward and reverse processes in Eq. 
become equal. These rates also will be large compared to 
the rates of collapse, transport, and composition change. 
The Ve fraction is defined to be Y^^ = {n^^ —Tip^)/ni,, i.e., 
the net number of electron neutrinos minus antineutrinos 
per baryon. The corresponding chemical potential is 

1/3 

fi^^ m 11.1 MeV(2/3ioyjy_, ) ' . In beta equilibrium we will 
have 

Me - A^i^e = A + Sninp, (2) 

where the difference of the neutron and proton kinetic 
chemical potentials is /i = fJ-n—fJ-p and the neutron-proton 
mass difference is 5m„p w 1.293 MeV. 

Neutrino trapping sets in for core densities pio ~ 10 
to 100, and the system attains beta equilibrium shortly 
thereafter. In the standard model of supernova physics, 
core collapse is adiabatic subsequent to neutrino trap- 
ping and the net electron lepton number (lepton frac- 
tion) Yl = Yy^ + Ye changes little. At neutrino trapping 
the electron fraction might be, for example, Yg « 0.35. 
However, as the collapse proceeds and the density rises, 
beta equilibrium shifts. By the time the core reaches nu- 
clear density and the collapse is halted ("core bounce"), 
the initial lepton number has been re-distributed between 
electrons and v^s, so that Yg ~ 0.30 and Y^^ « 0.05. 

The shock wave generated at core bounce will have 
an initial energy that is comparable to the in-fall kinetic 

energy and, therefore, scales as ~ Ye^^^ jsif. Note that 
any process which violates lepton number conservation 
and turns v^a into neutrinos of other flavors, including 
sterile flavors, will open phase space for electron capture. 
Electron capture prior to trapping and the establishment 
of beta equilibrium will lower Ye^ thereby lowering the 
initial shock energy. The additional electron capture on 
nuclei will increase the entropy to some extent. This 
will be off-set by post-trapping entropy loss if v,, — > Vs 
results in the production of sterile neutrinos which leave 
the core. 



B. The Neutrino Forward Scattering Potential and 
the MSW Resonance Condition 

An electron neutrino propagating through the stellar 
medium will experience a potential stemming from for- 
ward scattering on particles that carry weak charge (elec- 
trons, neutrons, protons, neutrinos). There will be no 
contribution to this potential from forward Ve ~ ^^^.^t scat- 
tering as long as the net muon and tau lepton numbers 
residing in the supernova neutrino seas are zero. Absent 
large scale neutrino flavor conversion, this is a good ap- 
proximation throughout the epochs of interest in the su- 
pernova environment. Electron capture produces z/g's, as 
outlined above, but i^e's and neutrinos and antineutrinos 



with mu and tau flavors can appear only through ther- 
mal pair production processes and these are suppressed 
by the large electron Fermi energies and low entropies 
characteristic of the in-fall epoch. Pair production of 
these species will be efficient after core bounce when the 
entropy is larger. 

With these provisos, we can express the z/g forward 
scattering potential (c/., Ref. (|j) as 

V = V2Gf (^Tie - in„^ + 2V2Gf (n^, - "pj 

+V2Gf {n., - npj + \/2Gf (71.. - npj (3) 
- + ,4) 

where Gp is the Fermi constant and = pNa is the 
baryon density, with p the density in gcm^'^ and where 
Njy is Avogadro's number and in the last equality we 
have set the net muon and tau neutrino number den- 
sities to zero. Here n^^ and rip^ are effective and Vf, 
number densities, respectively, so that n.^ — jt-p^ = nbY^^. 
Likewise, the net electron number density Ue and total 
(free plus nucleus) proton number density Up are related 
through overall charge neutrality, Up = Ue = n^- —11^+ = 
rtbYg. The total neutron density is then n„ = rib — Up. 

We will posit that in vacuum the unitary relation be- 
tween the propagating neutrino energy ("mass") eigen- 
states, \vi) and \v2), and the weak interaction (flavor) 
eigenstates is 

= cos9\vi) + s\i\0\v2) (5) 
\vs) = - sin 6'|i/i) + COS 6*11/2), (6) 

where \i'a) is the sterile neutrino state and 6 is an effective 
2x2 vacuum mixing angle. There would likely be mixing 
between the sterile neutrino and all three active neutrinos 
in vacuum, but in this work we will treat — Vs transfor- 
mations exclusively. We argue that this at least will serve 
to flesh out the main features of active-sterile conver- 
sion in collapse. This is probably a good approximation 
during in-fall because the small measured values of neu- 
trino mass-squared differences imply that active-active 
neutrino flavor conversion will be suppressed at densities 
characteristic of core collapse. Later, however, after core 
bounce, our treatment will be only approximate, as neu- 
trino flavor mixing in this regime will be 3 x 3 in nature. 
We can avoid the full 4x4 neutrino flavor mixing prob- 
lem because of the near maximal mixing between and 
Vr and the near identical interactions of these species and 
their antiparticles in the supernova environment pol l32| , 
at least until active-sterile and active-active neutrino fla- 
vor conversion has produced unequal and substantial net 
mu and tau lepton numbers. (See the discussion of the 
similar physics in the early universe in Ref. js^.) 

The effective in-medium 2x2 mixing angle 0m depends 
on the local neutrino forward scattering potential V and 
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can be found from 



sin^ 29m = 



A2 sin^ 261 



(7) 

(Acos26l-y)V A2sin2 26'' 

where, in terms of vacuum mass eigenvalue-squared dif- 
ference Sm^ = m2~mf and neutrino energy E^, we define 
A = 5m'^/2E^. The effective in-medium mixing angle is 
maximal (6*14 — 7i'/4) at an MSW resonance. A neutrino 
with energy Ei, and experiencing a forward scattering 
potential V will be MSW-resonant when 



A cos 20 = V. 



(8) 



The high densities encountered in stellar collapse dictate 
that at least one of the neutrinos must have a rest mass 
in the ~ keV range in order for there to be an MSW 
resonance for a typical (mean) neutrino energy |l4l |. As 
a result, we can safely set i5m^ « m| = m^, where nis is 
the rest mass associated with the sterile neutrino state. 
Furthermore, x-ray and closure constraints imply that 
the effective active-sterile mixing angle is small for this 
sterile neutrino mass range, so that we can set cos 26* « 1. 
With these approximations, the MSW resonance energy 
is 

« . ^„ f: — - (9) 



3V2GfpNa{Y, + Iy,^ - 
4.37MeV(ms/keV)^ 



P12 



±Y 

3 



(10) 



We follow convention and give the matter density p scaled 
by an appropriate value, so that pn = p/10"gcm~'^. 



C. Coherence and Adiabaticity 

Conversion of active neutrinos into sterile neutrinos 
can be facilitated either by active neutrino scattering- 
induced de-coherence or by coherent neutrino propaga- 
tion through an MSW resonance. The former process 
will dominate the production of sterile neutrinos at high 
density, where the scattering rate is large. For the neu- 
trino mass/mixing parameters of interest here, the latter 
process will dominate at lower density. Let us consider 
this process first. 

The fate of a neutrino propagating coherently through 
an MSW resonance depends on two factors. The first fac- 
tor is obviously coherence itself. By "coherence" we mean 
coherence on the scale of the MSW resonance width. 
Therefore, coherent neutrino evolution will correspond 
to large values (> 1) of the ratio of the neutrino's mean 
free path in medium to the MSW resonance width. 

The second issue in the neutrino's fate has to do with 
the degree to which neutrino flavor evolution is adiabatic. 
"Adiabaticity" is gauged by the ratio of the resonance 
width and the neutrino flavor oscillation length at reso- 
nance. Adiabatic evolution means that the neutrino re- 
mains in its particular superposition of initial instanta- 
neous mass eigenstate(s) as it propagates through reso- 
nance. 



Complete conversion of the neutrino's flavor will result 
if both of these ratios are large compared to unity and 
if the far asymptotic initial state of the neutrino corre- 
sponds closely to both an instantaneous mass eigenstate 
and a flavor eigenstate. For the densities and neutrino 
mass/mixing parameters of interest in stellar collapse, 
efficient coherent production of sterile neutrinos at an 
MSW resonance is always the result when the neutrino 
mean free path is large compared to the resonance width 
and the oscillation length at resonance is small compared 
to this width. 

The width of the resonance is St = H tan 26, where the 
effective weak charge scale height is H = \dlnV /dt\^^ , 
and where t can be either a time or length parameter 
along the neutrino's world line. The adiabaticity param- 
eter 7 for a neutrino with energy E^-bs is proportional to 
the ratio of St and the neutrino oscillation length at res- 
onance L°^g = inErcs/ {Sm'^ sin29), and is defined to be 
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27r- 



St 

T osc 



Sm'^n sin^ 29 



2Eic 



i29 



(11) 



Prior to neutrino trapping, v^s produced by electron 
capture will propagate at near light speed along radomly- 
directed trajectories. After trapping, the v^s will have 
nearly isotropic propagation directions. In either regime, 
however, sterile neutrinos i^s will freely stream at near 
light speed. The z^s's will be produced with randomly 
distributed propagation directions, reflecting the distri- 
bution of Ve propagation directions. 

Here we will consider sterile neutrinos which are di- 
rected radially outward. For this case, H computed along 
the radial direction is the relevant potential scale height. 
More tangentially-directed sterile neutrinos will experi- 
ence an effective potential scale height which is larger 
than that in the radial direction. As a result, our calcula- 
tions with radially-directed neutrinos will give lower lim- 
its on adiabaticity and, hence, will give under- estimates 
of neutrino flavor conversion efficiency. 

FigureHshows the potential scale height H for radially- 
directed neutrino trajectories as a function of radius. 
Also shown in this figure is the adiabaticity parameter 
7 relevant for sterile neutrinos with E^ — lOMeV and 
E^, = 100 MeV and with = 3 keV and = 10 keV, 
all for an assumed vacuum mixing angle in the Ve ^ i^s 



channel satisfying sin 29 



lo- 



in general, lower en- 



ergy neutrinos have smaller oscillation lengths at reso- 
nance and, consequently, larger adiabaticity parameters. 
This is evident in the figure. The "spike" in potential 
scale height and adiabaticity parameter evident at ra- 
dius r w 24 km in this figure, corresponds to a maximum 
in the potential V to be discussed below. 

After neutrino trapping, when v^s have relatively 
short mean free paths, the mean density experienced by 
a Ve will rise with time as the collapse proceeds. In this 
case, we can follow Ref. j3| ^-nd view the MSW resonant 
fiavor conversion in the mean rest frame of the Ve- In our 
one-zone collapse calculations we take the collapse rate 
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FIG. 1: The solid curve shows the eflective potential scale 
height Ti. (left vertical axis) in the radial direction as a func- 
tion of radius r. Also shown is the adiabaticity parame- 
ter 7 (right vertical axis) as a function of r for ^ Us 
with sin^ 26 — 10~^ and with sterile neutrino rest mass 
nis = 3 keV and energy Ei, = 10 MeV (long-dashed curve) 
and -Ei^ = 100 MeV (short-dashed curve), and for sterile neu- 
trinos with rest mass nis = 10 keV and energy E^, = 10 MeV 
(long dash-dot curve) and E^ = 100 MeV (short dash-dot 
curve) . 



to be, (see Appendix A) 

d In p , IS , 

^ = (loos-i)V;^, (12) 

and with this choice the adiabaticity parameter appro- 
priate for Ve — !■ Vs is 



At the very high densities encountered in the core near 
and after bounce, scattering-induced de-coherence can re- 
sult in significant production of sterile neutrinos (see for 
example the discussion in Ref. Q). The sterile neutrino 
emissivity in this case will be proportional to the product 
of the active neutrino scattering rate, a factor which takes 
account of the quantum Zeno effect, and sin^ 29m- This 
latter term will again make the sterile neutrino produc- 
tion rate largest at MSW resonances, where sin^ 20m = 1- 
The quantum Zeno factor is inversely proportional to the 
active neutrino scattering rate when this rate is large. 
This behavior tends to accentuate the enhancement of 
sterile neutrino production at MSW resonances. 



D. Evolution of the Potential 

We calculate the coupled histories of the neutrino for- 
ward scattering potential V, the electron fraction Yg, and 
the Ve fraction Y,y^ with a one-zone collapse/equation-of- 
state code described in Appendix A. This code includes 
a mean nucleus approximation and a self consistent cal- 
culation of entropy/temperature and electron capture 
rates on free nucleons and heavy nuclei. The density- 
time/radius histories of zones located at different radius 
will be self similar, so long as the collapse is homologous 
(in- fall velocity proportional to radius). With homology, 
all Lagrangian mass zones in the core will experience sim- 
ilar time-density histories, albeit with different beginning 
and ending points. In this case, the density of any zone 
at the instant of core bounce (which we take to be when 
the central density reaches pi2 = 3 x 10^) corresponds to 
the central density at an earlier time. Of course, electron 
capture and the associated fall in pressure will cause the 
actual collapse to deviate from homology. Nevertheless, 
our one-zone code serves to illustrate the qualitative de- 
tails important in understanding active-sterile neutrino 
flavor transformation and its effects on the equation of 
state in the core. 

The rise in density with time during collapse dic- 
tates the evolution of V and E^cs- The potential V 
rises and i?ros falls with increasing density {i.e., Eas. 1101 
and^. The MSW resonance in this case sweeps from 
higher toward lower energy through the distribution 
which, after neutrino trapping and the attainment of 
beta equilibrium, we can regard as Fermi-Dirac in char- 
acter, f,^ {E,) = [T^Fs (r/)]"'^^/ (g£„/T-„ ^ ^^ieK 
T is the temperature, 77 is the electron neutrino degen- 
eracy parameter p^^/T, and the relativistic Fermi in- 
tegrals of order k are defined by (77) = J-v+i ■ 
This is depicted in Fig. |5J If neutrino evolution is adia- 
batic, at any instant all the z^e's within the MSW width 
A£' — Ei-cs tan26' will be converted to sterile neutrinos. 

The one-zone calculation of Vf, forward scattering po- 
tential and resonance energy as a function of density for 
TUs = 3kev and sin^ 29 — 10^^ and for the constant col- 
lapse rate discussed above is shown in Figure |3I The 
striking feature in this plot is the peak in potential and 
corresponding minimum in resonance energy. 

The explanation for this peak/minimum is straightfor- 
ward. The trend of potential V with increasing density 
changes dramatically when the MSW resonance energy 
Ej-cs sweeps low enough to reach the Ve Fermi energy, fj,^^ . 
Given the degenerate conditions, there are not many v^^s 
with energies above the Fermi energy p^^. Early in the 
collapse, when E^^s > ^^vc^ not many sterile neutrinos 
are produced and there is only a modest decrease in Y^.. 
However, once i?res < Pv^i substantial numbers of v^s 
can be converted to sterile i/^'s and this can have a large 
effect on electron capture.. 

Neutrino flavor conversion i^e — i^s opens holes in the 
i/g distribution which allows for more electron capture 
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FIG. 2: The resonance energy Ercs sweeps from higher toward 
lower energy through the Ve distribution function (shown as 
a function of energy _E^^ ) , at any instant converting the 's 
contained in energy interval _Eres — Ai5 to iJrcs + Ai? (where 
AE = i?res tan 29) to sterile neutrinos Va. The number of i^e's 
converted at this instant is AA*'. 
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FIG. 3: One-zone calculations of Ue ^ Vs resonance energy 
Eras in MeV (left vertical axis, solid line) and forward scat- 
tering potential V in units of 10~*MeV (right vertical axis, 
dashed line) are shown as a function of density (in gcm~^). 
This calculation employs rris = 3kev and sin^ 29 — 10~^ 



and a concomitant decrease in electron lepton number 
per baryon Yi = Ye + Y-^^ . The portion of the change in 
lepton number stemming from the decrease in the num- 
ber of electrons, AKe, and the portion stemming from 
the decrease in the number of z^e's, AYi,^, is calculated 
self-consistently with the equation of state as described 
in Appendix A. Most of the lepton number decrease will 
be in Y^ , but some of it will be in Y^^ . In any case, the 
potential V is proportional to the factor Yg -t- "^v^ — \^ 
and this factor will decrease when there is substantial 
electron capture, e.g., when i?rcs < Mi^e- The net result is 
that V will decrease with increasing density beyond this 
point. 

In fact, there is a rough feed-back that develops once 
the resonance energy reaches the v^. Fermi energy. In- 
creasing density would tend to decrease i?res, but this is 
counteracted by a decrease in the factor Y^ -I- — \. 
The result is that the resonance energy tracks the 

1/3 

Fermi energy \x^j^ « 11.1 MeV(2pioyiyJ , increasing 
with increasing density and, more or less, staying just 
above /i^^. This is evident in Figure^ This figure shows 
-Eics as a function of density for a one-zone calculation 
with the same parameters as in Figure 13 but now also 
shows the corresponding values of the Ve Fermi energy 
[ly^ and the electron fraction Y^. 

The calculation of electron capture and, e.g., Y^ in the 
one-zone calculations depends on an accurate estimate 
of neutrino flavor conversion in the channel. 
In turn, as outlined above, this depends on the degree to 
which neutrino flavor evolution is coherent and adiabatic. 
There are three time scales in our one-zone calculation 
that give a rough idea of when these conditions obtain: 



Tcoib TLsc, and 8t 



res- ?coii is the mean-frec-time between 
I/qI^ is the oscillation time at resonance, 



collisions, Tqsc 
and <5ircs is the resonance width appropriate for core in 
fall as defined above. We then have 



TcoW — 



T 



f^n ^target ^rcl 



ml sin 26 ' 



5U 



■tan 26*. 



(14) 



(15) 



(16) 



Here (T„ is the cross section for a of energy E^^ to 
scatter on a neutron. We take this cross section to be 



1 



4 \meC 



where 



CTo 



4G2,m2ft2 



1.76 X 10-'''*cm2. 



(17) 



(18) 



In fact, the neutrino- nucleus coherent scattering cross 
section (crudely, larger than cr„ by nuclear mass-squared, 
A^) sets the mean free path, so our calculations here over- 
estimate Tcou- In Eq- H14|) . ritargct is the number density 
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FIG. 4: Same one-zone calculation as in Figure|3| resonance 
energy E^a& (solid line), electron fraction Ye (dashed line), and 
Fermi energy /i^^ (dotted line) are shown as functions of 
density. 

of nucleon/nucleus targets and Ui-d ~ c is the relative ve- 
locity between v^. and target nucleons/nuclei. Coherent, 
adiabatic flavor conversion Vp^ — + would require that 
the condition Tosc ^ ^^res <C TcoH be satisfied. 

The validity of this condition for the particular track 
of density, composition, and resonance energy in a one- 
zone calculation with ttIs = 3 keV and sin^ 26* = 10~^ 
is depicted in Figure [S] In this figure contours where 
Tosc/Tcou = 1, (5tres/Tcoii = 1, and T'coii/<5ires = 1 are 
shown as functions of resonance energy i^res and density 
p (in gcm~^). These ratios are < 1 below these contours. 
We conclude that coherence and adiabaticity is likely for 
all but the highest densities, near or beyond core bounce. 
Though these estimates are crude and leave out details 
of, e.g., detailed active neutrino cross section physics, we 
can conclude conservatively that neutrino flavor evolu- 
tion can be regarded as coherent and adiabatic for sterile 
neutrino parameters of interest for dark matter, at least 
through the epoch where the MSW resonance energy be- 
gins to track the Vf. Fermi energy. This epoch begins at 
the "knee" evident in the one-zone evolutionary track in 
Figure |5| This is the key issue in Kg evolution in the 
collapsing core. 



E. Electron Fraction Reduction 

The most striking change in core composition in the 
presence of i^e — ^ is the reduction in electron fraction 
Ye- In Figure El we show results of our one-zone calcu- 
lations with sin^ 29 = 10~^. This figure gives the final 
(at core bounce) value of Yg as a function of sterile neu- 




50 100 150 

Eres [MeV] 

FIG. 5: The one-zone calculation evolutionary track (solid 
line) with = 3keV and sin^ 26 = 10~^ is shown in the res- 
onance energy -Ercs, density p (gcm~^) plane. The conditions 
ToBc/Strcs = 1, Tosc/rcoii — 1 , and SUcs/TcoU = 1 obtain 
along the dashed, dotted, and dot-dashed contours, respec- 
tively. Values of these ratios are < 1 everywhere below these 
contours. 

trino rest mass (in MeV) for several different values of 
initial electron fraction y^""* and temperature T™'* (in 
MeV) at the onset of collapse. In broad brush, our cal- 
culations show that for sterile neutrino rest mass in the 
range 1 keV < nis < lOkeV and sin^ 26^ > 10"^^, a sig- 
nificant subset of the range of parameters of interest for 
dark matter, there is a substantial reduction in electron 
fraction, down to final values Yg « 0.22 from an initial 
value in the range y^'"'* = 0.35 to 0.4. This is in stark 
contrast to a one-zone calculation without Ve J^s, which 
would give a final electron fraction value Ye ~ 0.33 for 
the same range of initial electron fraction. 

Since, as outlined above, the homologous core at 
bounce forms the "piston" for the shock, the initial shock 
energy scales like ~ Ye^^^. As a result, the '-^ 30% reduc- 
tion in final electron fraction relative to the case without 
sterile neutrinos translates into more than a factor of two 
reduction in initial shock energy. 

And that is just the beginning of the problem. With a 
smaller homologous core there will be more material for 
the shock to traverse in the outer core. This is because 
the outer core is the remainder of the initial iron core 
which is not in the homologous, inner core. The entropy 
jump across the shock will be As ~ 10 units of Boltz- 
mann's constant per baryon, implying a shift in com- 
position in Nuclear Statistical Equilibrium from heavy 
nuclei in front of the shock to free nucleons and alpha 
particles behind it. This process is sometimes termed 
nuclear photo-dissociation. 

Since nucleons are bound in nuclei by ~ 8 MeV on av- 
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[MeV] 

FIG. 6: Final electron fraction Ye is shown as a function of 
sterile neutrino rest mass rUs as calculated in a one-zone code 
with sin^ 26 = 10~^ and initial values of electron fraction Fg""' 
and temperature T'"'' (in MeV) as shown. 



erage, the shock will lose 10 ergs per 0.1 Mq of outer 
core material traversed. Furthermore, our one-zone cal- 
culations show little electron capture-induced increase in 
core entropy during collapse relative to the no-sterile- 
neutrinos standard model. In part, this is because most 
of the extra electron capture stemming from holes in 
the Ve distribution opened by —^ Vs occurs when the 
medium is already in beta equilibrium. Also, sterile neu- 
trinos escaping directly from the core remove entropy. 
On balance, our one-zone calculations show only small 
net changes in entropy per baryon. 

The large sterile neutrino-induced decrease in initial 
shock energy and increased nuclear photo-dissociation 
burden may mean that shock stagnation occurs at a 
smaller radius than in standard core collapse models. In 
turn, with the shock stalled at lower radius, there would 
be less optical depth to neutrino absorption below the 
shock and, consequently, one would expect neutrino re- 
heating to be less efficacious. In short, the reduction in 
Yf, engendered by v,, — > Vs augers against a explosion. 

At this point one might be tempted to conclude that 
since (some) core collapse events lead to manifest explo- 
sions, the sterile neutrino mass/mixing parameters for 
which there is a significant Yg reduction are ruled out. 
This conclusion is ill advised for two reasons: (1) we do 
not know the detailed physics of shock propagation and 
electron capture/neutrino transport in the post-bounce 
regime; and, (2) sterile neutrinos can affect the rate at 
which lepton number and energy/entropy is transported 
to the region behind the shock. 



For example, based on a proper account of electron 
capture on very heavy neutron-rich nuclei in the post- 
bounce epoch, recently it has been shown that the stall 
point of the shock may be relatively insensitive to the 
value of Ye at bounce |3J| . This is despite the arguments 
to the contrary given above! The only robust conclusion 
that we can draw is that a large sterile neutrino-induced 
Yf, reduction will alter significantly the standard model 
of core collapse. These alterations would be on a scale 
far in excess of those stemming from nuclear physics and 
hydrodynamics/transport issues that modelers currently 
deal with. 

The second issue, sterile neutrino alteration of en- 
ergy/entropy/lepton number transport, is the subject of 
the next section. 



III. ENHANCEMENT OF ENERGY AND 
LEPTON NUMBER TRANSPORT 

The mechanism behind the explosion of core collapse 
supernovae is not well understood. The energy for the ex- 
plosion comes ultimately from gravitational energy, but 
the way this energy is transported seems to be crucial 
IMIMETi^, jalMEllillilEllil- contem- 
porary two and three-dimensional models show explo- 
sions E^ I41I Wh . while spherically symmetric, 
one-dimensional models, which boast the most sophisti- 
cated transport, equation of state, and nuclear physics 
do not seem to yield convincing explosions. Perhaps this 
reflects the actual situation in nature. However, many 
uncertainties in the fundamental physics input remain in 
these models and it is an open question whether vigorous 
explosions will occur in one dimension [49l |. 

Simulations suggest that a modest increase in neu- 
trino energy luminosity from the proto-neutron star sur- 
face (neutrino sphere) might increase neutrino energy 
deposition rates behind the shock enough to cause a 
viable explosion in what would otherwise have been a 
dud. This shows that neutrino energy transport in the 
proto-neutron star core may be the crux issue in the one- 
dimcnsional core collapse explosion mechanism problem. 

The peak in the potential V and associated dip in the 
resonance energy E^-es found in our calculations suggests 
a process whereby the Vf, neutrino luminosity might be 
increased. With the local minimum in E-^ca near the mid- 
dle of the core, a neutrino of a given energy could en- 
counter two MSW resonances: one deep in the core and 
one further out, closer to the neutrino sphere. Figure [7| 
shows the resonance energy and radius parameter in the 
in-falling, pre-bounce core as functions of density. 

Imagine that there is a source of sterile neutrinos from 
flavor conversion occurring deep in the core, where v^. en- 
ergies are ~ and so are large, perhaps ~ 100 MeV. 
This sterile neutrino emissivity could arise from medium- 
enhanced coherent flavor conversion at the in- 
ner MSW resonance or, more likely, scattering-dominated 
de-cohrerent flavor conversion in any channel Va t^s 
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FIG. 7: Resonance energy (left) and radius parameter r 
(right) in the in-falling, pre-bounce core are shown back-to- 
back as functions of density p (vertical axes). An example 
given resonance energy E-^cs corresponds to two locations, ri 
and r2, and two corresponding densities, p\ and p2. 

{a = e, /i, r) or medium-enhanced, de-coherent v^, 
occurring at the inner MSW resonance. 



The shock is generated at the edge of the inner core and 
will alter the entropy and, hence, the composition in the 
outer core. 

However, some gross features of our shock energy en- 
hancement scheme are easy to discern. If, for example, 
the feedback mechanism relating z/e Vs flavor con- 
version and potential V continues to hold through the 
epoch when the shock propagates through the outer core, 
then, as in the in-fall epoch, the resonance energy E-^cs 
will track the chemical potential /i^^ and will remain 
slightly above it, i?res ^ Mi^c- If in addition, /i^^ sim- 
ply tracks the density so that near the center of the core 
~ 150 MeV, then nearly all the sterile neutrinos cre- 
ated by scattering-induced de-coherence deep in the core 
will be re-converted to active v,, neutrinos further out. 
Taking the high energy Vs emissivity from Ref . and as- 
suming complete conversion of this energy into fe's which 
thermalize quickly nearer the core's edge, we estimate 
that < 10^^ ergs of additional neutrino energy could be 
radiated from the neutrino sphere during the ~ 1 s dura- 
tion of shock re-heating. This is, however, a very crude 
estimate. We will estimate the effects of shock propaga- 
tion on the neutrino flavor transformation potential in a 
subsequent work. 




FIG. 8: High energy VeS could be converted to sterile neutri- 
nos deep in the core and then re-generated as further out, 
nearer the neutrino sphere (edge of core). 

These high energy sterile neutrinos would move out at 
near light speed and encounter the second, outer MSW 
resonance. There they could be coherently re-converted 
Vs Vf, to active VeS. The result would be an upward 
re-normalization of the net electron lepton number and 
energy transport rates. In other words, flavor mixing 
would allow u neutrinos to spend part of the time in 
sterile states that have higher effective diffusion rates. 
This is depicted in Figure |S1 

For this mechanism to aid the shock re- 
heating/explosion process , the rough arrangement 
of an inner source of sterile neutrinos at high energy 
and an outer MSW resonance would have to survive 
shock propagation through the core. This is because 
the epoch at which shock re-heating occurs is between 
a few hundred milliseconds and 1 s after core bounce. 



IV. DISCUSSION AND CONCLUSIONS 

We have demonstrated that the existence of sterile 
neutrinos with rest masses in the range of interest for 
dark matter could have significant implications for the 
in-fall, collapse epoch of core collapse supernovae. In 
particular, sterile neutrinos with rest masses in the range 
IkeV < rUs < lOkeV, and possessing effective vacuum 
mixing angles characteristic of the channel Ve ^ Vs which 
satisfy sin^ 261 > 10~^^, could result in enhanced elec- 
tron capture during in-fall with a concomitant decrease 
in electron fraction Yg- As discussed in Section II E, the 
decrease in Y^. will translate into a decrease in relativis- 
tic electron degeneracy pressure in the core. In turn, the 
consequence of this will be a smaller homologous core 
and, therefore, a lower initial shock energy and an in- 
creased photo-disintegration burden on the shock. Al- 
though these alterations certainly auger against obtain- 
ing a viable shock and explosion in the framework of cur- 
rent models, it must be kept in mind that there remain 
significant uncertainties in fundamental nuclear electron 
capture, equation of state, and transport physics in the 
post core bounce regime. 

However, we could conservatively conclude that the 
existence of sterile neutrinos with these properties would 
result in a supernova paradigm radically different from 
current models. Additionally, we have outlined a means 
whereby active-sterile-active neutrino flavor transforma- 
tion (specifically, v^, —^ — > Ve) could alter electron 
lepton number and energy transport rates. Active neu- 
trinos which spend part of their time propagating as ster- 
ile neutrinos will have larger transport mean free paths 
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FIG. 9: Lightly shaded and heavily shaded regions show the 
ranges of sterile neutrino rest mass nis (in keV) and effec- 
tive Ue ^ i^s vacuum mixing sin^ 26 that give large Ye re- 
duction and shock re- heating/explosion enhancement, respec- 
tively. The line labeled by 7100 > 1 shows the contour of 
mixing parameters along which the adiabaticity parameter 7 
at the outer MSW resonance is unity or greater for neutrinos 
with energy Ei, = 100 MeV. The corresponding line 710 > 1 
shows where 7 is unity or greater for neutrinos with energy 
-E^ = 10 MeV. The heavy solid line and the cross-hatched 
region to the right show the range of parameters ruled out by 
X-ray observations of the Virgo cluster 0, . The dotted 
line and the region to its right is ruled out by observations of 
the diffuse X-ray background 50]. The dashed line (labeled 
LMC) is the claimed X-ray limit from the Large Magellanic 
Cloud [sH. The dashed and solid horizontal lines are the 
Lyman-a forest-derived lower limits on iris from Ref. and 
Ref. IH^], respectively. 

on average. In short, much of supernova physics "down 
stream" of core bounce could be quite different from stan- 
dard models. 

For example, a higher i/g luminosity as a result of 
j'e — > I's — > i^e would result in a number of changes in 
the near hydrostatic envelope above the neutrino sphere 
and below the shock. This would mean a much reduced 
neutron excess in this region, since the neutron-to-proton 
ratio is determined locally by a competition between 
i^e+n ^ p+e~ and Pe+P ^ n+e'^ . This would have im- 
plications for alpha-rich freeze-out, neutrino-heated nu- 
cleosynthesis coming from this region. In fact, this effect 
could be beneficial, increasing Yg and lowering the yield 
of neutron-rich, neutron number TV = 50 nuclei. How- 
ever, the initial neutrino spectra at the neutrino sphere 
could be altered by the the enhanced deposition of Ve's 
just below the proto-neutron star surface. This could 
partially counteract the neutron excess-lowering effect in 
the envelope by reducing the average Ve energy relative 
to that for P,,. 



In the end, our results call into question the old idea 
that sterile neutrinos remove energy from the core and, 
therefore, cannot exist, else there would be no super- 
nova explosions. It must be kept in mind that the super- 
nova explosion energy (kinetic plus optical) is ~ 10*^^ erg, 
which is only ~ 1% of the total amount of energy radi- 
ated away over some 10 s as neutrinos of all kinds. It may 
be possible to "throw away" much of this energy and still 
get an explosion. The timing, rate, and location of en- 
ergy deposition are the most important determinants of 
explosion physics and we have demonstrated that some 
or all of these can be affected by sterile neutrinos with 
masses/mixings in the range of interest for dark matter. 

In Figure El we summarize the the parameter space of 
sterile neutrino rest masses and effective vacuum mixing 
angles (in the z/g ^ channel) that gives rise to super- 
nova effects. This figure also shows possible constraints 
stemming from cosmological sterile neutrino radiative de- 
cay — > t'ct + 7) . To obtain significant Ve i^s i^e 
enhancement of luminosity at the neutron star surface 
(neutrino sphere), high energy z/^'s coming from deep in 
the core must encounter the outer MSW resonance and 
be converted with high efficiency to i^e's. In the figure, 
the line labeled by 7100 > 1 shows the contour of mix- 
ing parameters along which the adiabaticity parameter 7 
at the outer MSW resonance is unity for neutrinos with 
energy = 100 MeV. Likewise, this figure also shows 
the corresponding contour, 710 > 1, for adiabaticity for 
neutrinos with energy Ei, — 10 MeV. The lightly shaded 
region in this figure, which shows the parameter range for 
large reduction in Ye , is bounded on the low mixing angle 
side by the contour 710 > 1. This is because augmenta- 
tion of electron capture requires only that neutrinos 
with energies E,y ^ 10 MeV be efficiently converted to 
j/s's. 

The existence of relatively small scale structure in the 
universe, i.e., the Lyman alpha forest, sets an upper limit 
on neutrino suppression of the matter power spectrum 
and the dark matter particle coUisionless damping scale. 
In turn, this sets a lower limit on the rest mass of sterile 
neutrinos if they are to be the bulk of the dark mat- 
ter. This limit is in dispute. In Figure]^ we show both 
the Abazajian limit |53i | and a more conservative limit 
[s^ l. Likewise, the lack of an observed X-ray line feature 
in observations of the Large Magellanic Cloud (LMC) 
has been argued to set an upper limit on sterile neu- 
trino rest mass '5l|. We also show current diffuse X-ray 
background constraints jH^]- Perhaps the most stringent 
constraints come from the lack of observation of a sterile 
neutrino X-ray decay line from the Virgo cluster 4, 2§l. 
This is the heavy solid line and the cross-hatched region 
to its right. 

With all constraints applied, there remains a signifi- 
cant region of unconstrained TOs— sin^ 29 parameter space 
in which there could be large supernova in-fall epoch and 
shock re-heating effects. Neutrinos with these parameters 
could also be the dark matter, at least if there is a signif- 
icant primordial lepton number Q, ^ |^ or with par- 
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ticular couplings |2] ■ These neutrinos fall in the range of 
parameters that, after core bounce, could engineer large 
pulsar kicks pjl Il2| | and modify hydrodynamic neutrino 
energy transport to the base of the shock • 

The ultimate conclusion that can be drawn from our 
work is that stellar collapse is remarkably sensitive to new 
lepton number- violating physics in the neutrino sector. 
The reason for this is clear: the energy in core collapse 
supernovae is gravitational. This gravitational binding 
energy is stored in degenerate seas of leptons {ve's and 
electrons). Any lepton number violating processes in op- 
eration can tap into this energy and thereby alter super- 
nova core composition and dynamics. Sterile neutrinos 
with rest masses ~ keV which mix with active neutri- 
nos at the level of one part in ^ 10^^ in vacuum are 
likely to be undetectable in conventional terrestrial ex- 
periments but, remarkably, could completely change the 
current model for core collapse supernovae. 



The mean nuclear mass is A. This is determined by a 
minimization of nuclear energy which gives the condition 
Wsurf^^^^ — 2Wcouiomb^^^^, from which we obtain 

A = 194(1 - Y,f{l - 0.234pif -I- 0.00194pi2)"^ (A5) 

With these expressions it is evident that /U„ and /t are 
function s of Yg , Xn , and p 

fi^fiiY,,Xn,p), (A6) 



Pn = Pn{Ye,Xn, p). ( A7) 

Likewise, as discussed in section II, the electron and 
chemical potentials are functions of p, temperature T (or 
in energy units fcsT, where ks — 0.08617 MeV/Tg is 
Boltzmann's constant), and Ye or Y^y^, 

Ml/, = /ii/e (Y^. 7 ksT, p), (A8) 
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APPENDIX A: NUMERICAL CALCULATIONS 

To gauge the effects of altered electron capture physics 
in the in-fall epoch of stellar collapse, we employ a mod- 
ified version of the one-zone code described in Ref. [3ll |. 
For the neutron kinetic chemical potential (in MeV) we 
use 



fin = -16+125(0. 5-r,)-125(0.5-r£)^ 



Wsurf i-7Y 



^1/3 2(1 - Y,) ' 
(Al) 

while for the difference of the neutron and proton kinetic 
chemical potentials, ft = Hn — l^p, (again, in MeV) we 
employ 



A* 



250(0.5 - Y,) - W.^.iA-^/'' ( + 



2 1-21; 



Y, Y, l-Y, 



(A2) 

In these expressions, Y^ = ^e/(l — A„), where A„ is 
the free neutron mass fraction. Following Ref. ^3^, the 
nuclear surface energy Wsmf and the nuclear Coulomb 
energy Wcouiomb are taken to be: 



Wsurf = (290MeV)y,2(i_y^)2 



(A3) 



W^couiomb = (0.75MeV)y,2(i _ 0.234^: 



1/2 



0.00194/912). 

(A4) 



p.e = p.e{Ye, kBT,p). 



(A9) 



The chemical potentials of the leptons and the free nu- 
cleons are related through the beta equilibrium condition, 
Eq. (O . We assume that beta equilibrium is attained on a 
time scale short compared to the dynamical time implied 
in Eq. (|12(l . The beta equilibrium condition is then a 
constraint among Ye, Y„^, ksT, and X„, which we could 
describe as 



//3-oqil(A„, ksT, Ye,Y^^,p) — 0. 



(AlO) 



Note that we take the time evolution of p explicitly from 
Eq. 

Another constraint on the system can be obtained 
by assuming Nuclear Statistical Equilibrium (NSE), in 
which case the Saha equation gives 



A„ 



™(Ml!!!„p(^) 

Pw \kBi J 



(All) 



This is a good approximation whenever nucleons are non- 
degenerate, i.e., well below nuclear saturation density at 
pi2 ~ 300. The assumption of NSE provides a constraint 
among Ye, ksT, and A„, and could be written. 



/nse(A„,A;bT, Ye,p) = 0. 



(A12) 



Since particle number in both election capture and 
active-sterile flavor transformation Ve Vs is preserved, 
we have an additional constraint on the number of spin- 
1/2 leptons. 



Ye + Y„. 



Y„ 



Yl 



const, 



(A13) 



where Y^ is the initial (at neutrino trapping) lepton num- 
ber per baryon. Equivalently, we can consider the time 
derivative of this constraint equation. 



Ye + Y,, 



Y,, 



0. 



(A14) 
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The sterile neutrinos are not trapped, of course, so Fj^, 
could be thought of as a sterile neutrino number emissiv- 
ity per baryon. 

The entropy per baryon in our calculation is regarded 
to be a funcion of Yg, ^VaJ -^n, ksT, and p, namely 
S{p, Ye, Y^^ , Xn, ksT). Therefore, the entropy evolves in 
time according to 

= ^(pW,yi"\yJ:\xi"\fcBT(")) + SM, (A15) 

where At is a time increment related by Eq. (|12|l to the 
increment in density. We could write the total differential 
of entropy in a time/density step as 



dS 



1 



- (^flgdYe + (pn + mnC^)dYn 

+ {jjip + mpC^)dYp + ^v^dYi,^ 
1 



ksT 

-(/ie - - /i 

1 



)dYe 



knT 



[dQ + dYi,^ ] 



(A16) 



(A17) 
(A18) 



where dQ^^^^^ and dY^^ are the heat loss per baryon and 
stelile neutrino number generation, respectively, stem- 
ming from the flavor conversion — > i^s occuring during 
the time step as the resonance energy i^ros (Eq. ((117)1 1 
sweeps through a portion of the Ve distribution. Note 
that the /3-equilibrium condition (imposed to simplify the 
second-to-last equality above) gives a simple relation be- 
tween the time rates of change of the entropy and heat 
loss per baryon and the sterile neutrino production rate 
per baryon. 



5" = 



1 



knT 



(A19) 



At an MSW resonance with physical width St^^s and total 
Ve neutrino energy contained in this width l^E^^ {Yi^^ , T) , 
we have 



AE,SY,^,T) 



,{Y.,,T,p), (A20) 



Combining the constraints in Eq. I|A13|) and Eq. (|A21|) 

we obtain 

yjn+l)^yj„+l) =yJn)+yW_l> (FJ") , /cfiT^") , p^") ) At. 

(A23) 

The other two constraints, Eqs. lAlOl and IA12I can be 
expressed as 

W(x("+^\4'^+'^fcsT("+i),i;("+i),yj:+i),p("+i)) = 0, 

(A24) 

and 

/NSE(X<i"+'\fcsT("+i),yi"+i),p("+i)) = 0. (A25) 

Eqs. (|X22|l . (|X23)l . l|X24)l . and l|J25)l form a set of non- 
linear equations with respect to variables fJ"^^', yJ)^~^^\ 
and kBT'^''+^^ and can be solved numerically, 
e.g., by a Newton-Raphson method. 



= = (.Y,. ^T,p), (A21) 

otrcs(P) 

where AA'',^^ (Yg, Y^^ , T) is the number of v^s within the 
resonance width. With these constraints, definitions, and 
simplifications, Eq. l|A15p becomes 

5(^(n+l) ^ yjn+l) ^ yjn+l) ^ ^(n+l) ^ (n+l) ) 

= 5(p("),ri"),yH,x("),A:BT(")) 

+S{yI;^^ , keT^''^ , ) At. (A22) 
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